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Abstract 

The non-singlet and singlet anomalous dimensions of the twist-2 light-ray operators for 
unpolarized and polarized deep inelastic scattering are calculated in 0(a s ). We apply these 
results for the derivation of evolution equations for partition functions, structure functions, 
and wave functions which are defined as Fourier transforms of the matrix elements of the 
light-ray operators. Special cases are the Altarelli-Parisi and Brodsky-Lepage kernels. 
Finally we extend Radyushkin's solution from the non-singlet to the singlet case. 
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1 Introduction 



The study of the Compton amplitude for scattering a virtual photon off a hadron is one of the 
basic tools in QCD to understand the short-distance behavior of the theory. In the general case 
of non-forward scattering the Compton amplitude is given by 

T, v (p + ,p.,Q) = % J d 4 xe^(p 2 \T(J,(x/2)J u (-x/2))\ Pl ), (1) 

where p± = p 2 ±pi, Q = (qi + qz)/2, and p% + q 1 = p 2 + q 2 . 

The time-ordered product in eq. ([!]) can be represented in terms of the non-local operator 
product expansion [jj : 

T{J lx {x/2)J v {—x/2)) « / ck_ / dn + C a (x 2 , k_, k+, fi 2 )S^x p O"(K + x, k_x, /j 2 ) 



OO J — oo 



+ C a ^(x 2 , K+, fi 2 )e flu pa XpO^ (T (K + X, K_X, /i 2 ) (2) 

with S^xt = g^gxc - fi^AfiW - g^gxv and e pupa denoting the Levi-Civita symbol. The light-like 
vector 

x = x - r](x.r] T].r]) + r] Jl - — (3) 

77.77 y (x.ri) z 

is related to x and a constant subsidiary four-vector 77 which drops out in leading order. Un- 
like the local operator product expansion 0, which is usually applied in the case of forward 
scattering, eq. (fj) straightforwardly leads to compact expressions for the coefficient functions 
C a (x 2 , k±, fi 2 ) and light-ray operators 0^(k±x, /i 2 ). Indeed this expansion is a summed-up lo- 
cal light-cone expansion [|l], In lowest order it contains only quark operators with two external 
legs. These operators are to be decomposed into operators of different twist. Here, we consider 
the contribution of twist-2 operators onlyQ for which we calculate the anomalous dimensions 
entering the corresponding renormalization group equation |J. 

Different processes, such as deep inelastic (forward) scattering ||, deeply virtual Compton 
scattering [5,7-12], or more generally non-forward scattering processes, cf. e. g. JI3J, and pho- 
toproduction of mesons [13-16] contain as an essential input the same light-ray operators of 
twist-2: All interesting non-perturbative partition functions are Fourier transforms of corre- 
sponding matrix elements of these operators. Therefore, all evolution equations result from the 
renormalization group equation of the light-ray operators involved. The evolution kernels can 
be determined from the anomalous dimensions of these operators. This procedure will be car- 
ried out here. The Altarelli-Parisi kernels as well as the Brodsky-Lepage/Efremov-Radyushkin 
kernels appear as special cases. 

A second interesting feature of this approach is that non-forward processes contain more 
variables than forward processes. In the generalized Bjorken region this implies that instead of 
only one distribution parameter in the forward case, the Bjorken variable x, two distribution 
parameters appear in the present case. One of them is related to the sum and the other to 
the difference of the external momenta. The explicit solution of the evolution equation in the 
non-singlet case was firstly found by Radyushkin jL7| . Here, we extend his method to the singlet 
case. On the other hand, by imposing specific conditions on the momenta also one-parameter 
equations may be obtained which are equivalent to the general case treating a parameter which 
is fixed for special free second variable. 

2 The notion of twist is not quite unambiguous in this context. It is rather used to label the type of operators 
being concerned, for which in the limit P2 — > pi the notion applies Mj. 
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In this paper we present the results of a calculation of the non-singlet and singlet twist- 
2 anomalous dimensions for the general case both for unpolarized and polarized scattering, 
cf. 0. We derive general evolution equations both for the operators and matrix elements in the 
non-forward case. Furthermore we derive a series of special cases previously discussed in the 
literature. 



2 Twist-2 Light-Ray Operators and Different Choices of 
Partition Functions 

We consider the following twist-2 flavor singlet operators^ which are known to mix under renor- 
malization 

2 r -| 

9 (ki,K 2 ) = - ^(KiX^X^r^x) - Vv(/S2£)7^Vv(/Sl£)J , (4) 

G ( Kl , K2 ) = x^'F a l{K X x)F\ lv {K 2 x), (5) 

1 r -i 

Ol{Ki,K 2 ) = - [^r(KlX)>y 5 'y ll xV'ljj r (K 2 x) + V , r(«25)757/i £A 'V'r(«l£)J , (6) 

0%(k u k 2 ) = ±&>x*' [^ a ;( Kl j)FV',(«2x) - FV(«2^)^V i '( K i i )] • (7) 

Here ip r denotes the quark field, F" and F^ v the gluon field strength and its dual, respectively, 
and furthermore we use 

K\ = K+ — k_, k 2 = k + + k_ . (8) 

To simplify the considerations we apply the axial gauge, x^A^ = 0, by which the phase factor 
connecting gauge dependent fields equals to unity. 

The operators ([5]) and (0) appear only in the next order in a s in the Compton scattering am- 
plitude. However, due to the mixing of operators of the same dimension and twist they emerge 
in the renormalization group equations. Whereas for the determination of the anomalous dimen- 
sions it is sufficient to investigate the operators (§-0), by considering the evolution equations 
for observables it is necessary to introduce generalized partition functions related to the matrix 
elements of the operators O l ,i = q,G; (the index 5 will be suppressed in the following). The 
expectation values may be represented by either 



(ixp + ) hi 

or 

r+l 



dz + dz_e- iK - {ip+z++ip - z - ) F i (z+, zJ) (9) 



(Pi |0*|p2> • (K-) hi = e~ iK+Sp - J dz + dz.e- iK -^ p+z++ ^- z - ) G i {z + ,z^) . (10) 

Here hi denotes the degree of homogeneity, hi, with respect to a rescaling of the k- variables of 
the quarkonic and gluonic operators (cf. [[jj) with 

h q = l, h G = 2. (11) 



3 In leading order the non-singlct operator is related to eq. (Q). The corresponding results which are derived 
below for this operator apply literally to Ons(ki, K2) in this order. 
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One way to express the partition function, eq. (§), is connected with the states and has the 
property, that in the forward limit, \pi) — > |p 2 ) , these amplitudes lead to the well known structure 
functions appearing in the Altarelli-Parisi-equations. The second expression, eq. (p70|) , which is 
related to the light-like spread of the operators leads to simpler two-variable evolution equations 
the solution of which is more straightforwardly obtained, cf. section 7. 



3 Anomalous Dimensions of Twist 2 Operators 

The renormalization group equations for the singlet operators read: 



d/j 2 



2tt 



DaK l3 {a\ : a 2 , k^Oj^, k' 2 ) 



with 



a s 



Da = da±da 2 0{1 — «i — a 2 ) #( a i) 0(a 2 ), 
g 2 /(An) the strong coupling constant, \i the renormalization scale, and 



:i2i 



(13) 



k[ = (1 — + ttiK 2 , n' 2 = a 2 Ki + (1 — ol 2 )k 2 . 

By K = (K^) we denote the matrix of singlet anomalous dimensions in the unpolarized case. 
The singlet evolution equations for the polarized case are obtained by replacing the operators 
qq,g ky qq,g an( j j-^g m atrix K by AK. As far as relations are concerned which are valid both 
for the unpolarized and the polarized case under this replacement, we will, for brevity, only give 
that for the unpolarized case in the following. The matrices of the singlet anomalous dimensions 
are 



K 



K gg K iG 

K Gq K GG 



and AK 



AK qq AK qG 
AK Gq AK GG 



(14) 



respectively. In leading order, the non-singlet anomalous dimensions obey f^ NS = AK^ S 
K qq = AK qq . 

For the unpolarized case the anomalous dimensions are given by 



K qq (a h a 2 ) = C F y. - 8(on) - 6(a 2 ) + 6(ai) 

K qG (a u a 2 ) 
K Gq ( ai ,a 2 ) 
«i, a 2 ) 



1 

a 2 

— 2A^/Trk_ {1 — a\ — a 2 + 4aia 2 } , 
-C F —{5( ai )5(a 2 ) + 2}, 



+ S(a 2 ) 



Ct\ 



-5(ai)5(a 2 ) 



K GG 



C A {4(1 — oc\ — a 2 ) + 12aia 2 
1 



(15) 

(16) 
(17) 
(18) 



+ 5(a 1 



ct 2 



2 + a 2 +S(a 2 ) 



a>i 



-2 + ai 



+ y^(ai)^(«2) 



in 0{a s ) where C F = {N 2 C - 1)/2N C = 4/3, T R = 1/2, C A = N c = 3, /3 = {11C A - 4T R N f )/3. 
The [ ] + -prescription is defined by 



dx [f(x,y)],ip(x) 



dxf(x,y) [ip(x) - if 



(19) 



if the singularity of / is of the type ~ l/(x — y). 
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Correspondingly the anomalous dimensions for the polarized case are 



Aif«(ai,a 2 ) = K^( ai ,a 2 ), (20) 
AK qG (a h a 2 ) = -2N f T R K- {1 - a x - a 2 } , (21) 

AK Gq ( ai ,a 2 ) = -C F —{5( ai )6(a 2 )-2}, (22) 

AK GG ( ai ,a 2 ) = K GG (a ll a 2 )-12C A qna 2 . (23) 

Whereas the anomalous dimensions for the polarized case are derived for the first time in 
those for the unpolarized case were found several years ago [0, §]. Recent calculations have also 
been performed in || |1(| . The matrices K % i and AK l i determine the evolution of the operators 
Of£,Ol 5) , and 0? 5 \, respectively, in 0[ 



O'.s 



The determination of these anomalous dimensions is straightforward. In lowest order one has 
to calculate the one-particle irreducible one-loop Feynman diagrams containing the operators 
considered as first vertex. The corresponding Feynman rules and the calculation in the covariant 
gauge for one example are presented in ref. J7|. Here we have performed the calculation in 
the axial gauge which leads to essential simplifications. We used the Leibbrandt-Mandelstam 



prescription for the axial poles [IS . 



4 Evolution Equations of Partition Functions 

The non-perturbative partition functions defined in section 2 have to be determined experimen- 
tally at an input scale Qq. Their values at higher scales Q 2 can be calculated perturbatively 
solving the corresponding evolution equations. Here we derive only the equations which emerge 
from our two choices eqs. (|9|Jl0f). These equations follow immediately from the renormalization 
group equations of the operators, eq. fljjl). 

For simplicity we introduce different variables which are most convenient, cf. |7j], for the 
subsequent discussion 

K + ~ K + i K ~ /CA\ 

Wi = cti — a 2 = , w 2 = 1 — ai ~ ol 2 = , (24) 

Da = \Dw = ^dw\dw 2 Q(l — w\ — w 2 )Q{l + w\ — w 2 )Q{w 2 ). 
In these variables the RG-equation ( |l~2"|) reads 

fx 2 -^0\K + , «_) = ^4^- [ DwK»{ Wll w 2 , K^)O j (n + + twi«_, tw 2 «_). (25) 
djji z 2n J 

By Fourier transform we get an evolution equation for the distribution functions G x , inserting 
the definition ([R]) also into the right hand side 

2_^_ r i( \_ [ + °° d.K-Xp+ r+°° dK^0Cp_ i K _(xp-Z-+xp + z+) ( \hi 



dfi z J-oo 2n J-oo 2n 



An 



J DwK ij ( Wl ,w 2 ,K_) J J dzidz , + (Ki)- h ^e- iK '- {ip+z '+ + ^- z '- ) G j (z' + z'_). 



4 A different, but related, quantity was studied in ui 
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The ^--dependence of the evolution kernel is K^(wi, w 2 , t_) = K^(wi, w 2 ){K>-) hi hi . Carrying 
out the integrations with respect to dK_xp± we get 

8(z- — W2z'_ — wi)5(z+ — w 2 z' + )G J '(z' + z'_) (26) 

leading to 

^G'(z + ,z_) = «f / +1 r 1 fefoV+VJ . (27) 

U/i 47T J — 1 J — 1 i-^+i ^+ ^+ v ^+ / 

The derivation of the evolution equations for the partition functions @ follows the same line. 
The difference is that the integration with respect to dn-xp + leads not directly to a 5-function 
so that one obtains 

^-^- 2 F\z + ,zJ) = £Y_r dz'_dz' + r^(z + ,z^z' + ,zL)F^z' + ,z'_), (28) 

where 

P»(z+,z_,<.,0 = /= / d J^P± e _ K^z_-w 2 z'_,w 2 ) , (29) 

with /ijj = hi — hj. Wether the <i/€_xp + -integration leads to a 5-function, its derivative, or an 
integration depends on the values of i and j. In the case of an integration a 9-function appears, 
for which the integration constant has to be determined separately. 



5 One- Variable Evolution Equations 



The equations given above cover the most general case depending on two partition variables z + 
and Z-. To reveal the relation of the general case to more special cases we introduce a kinematic 
condition by xp^ = rxp + . This condition appears natural having in mind scale invariance. For 
special cases being discussed below the parameter r will be fixed. On the other hand, one may 
consider r as a second general parameter aside xp + resembling the general two-variable case. 
Instead of the variables z_ and z + we use t and r as new variables defining 



Fi(z + ,z_) 



/+oo 
dz_Fi(t — tz_, Z-), 
-00 

I f+oo f+oo \z \ i \ , 

dt dT L^ e ™-{* + -t)/T^ T) 

{2tI) Z J-oo J-oo t z 



(30) 
(31) 



These partition functions are related to the expectation values of the operators Q q and 0° by 

{pi\O q (-K^X, K_x)\p 2 ) 

XP— =TXpj r 



(ixp+) 

(pi\0 G (-K-X, K-X)\p2) 



-CO 



(ixp + y 



-oo 
+00 



dte- iK - ip+t <5> q (t,T), 



(32) 
(33) 



Xp—=TXp+ 
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The evolution equations are : 



dt'v:ut,t',T)&(t',T). 



dp, 2 y ' ' 2tt 
The corresponding extended kernels read 

" f 1 f + °° _ 

Vext(t,t',r) = J DaK lJ (a 1 ,a 2 ,K-)— J ^ d(p + xn^) [ip+x] lJ 

x — exp {ij9 + /t_a; [i — (1 — «i — a^H' + t(«i — o^)]} 
and obey the scaling relation 

For convenience we write the general expressions for the evolution kernels in the variables 



x 



t r 



They are given by : 

Vext{ti t j r 

with 



^ { V™(s, y) + y)p(x, y) + \C F 8{x - y) } 

^ {^ G (:r,y)p(:r,y) - ^ G (x, y)} f2y " r 
^ { \/ G «(x, y)p(x, y) - F^x, y)p(x, y)} 
^- { V GG (:r, y)p(x, y) + V GG (x, y)p(x, y)} 



2 
2 

2a; - 1 
2y-l 
2a;- 1 



+—(3 5(x - y) 



V e q x q t (t,t',r) 

± {AV> G (x, y)p(x, y) - AV G (x, y)p(x, y) } 
^ {AV°*(x, y)p(x, y) - AF G9 (x, y)p(x, y) } 
{ AV GG (x, y)p(x, y) + AV GG (x, y)p(x, y) } 



2r 
1 

27 



2y-l 
2 
2 

2x - 1 
2y-l 
2a;- 1 



+—P 5(x-y), 



p{x,y) = e^.-^e^j sign(y) 



and 



V qq (x,y) 



Ci 



x 1 
+ 



y y (y - x) 



V* G (x,y) 
V Gq (x,y) 

V GG (x,y) 

AV qq (x,y) 
AV qG (x,y) 

AV Gq (x,y) 
AV GG (x,y) 




y + x 



(48) 
(49) 

(50) 

(51) 
(52) 

(53) 

(54) 



Note that the kernels given in eqs. (j38|- |5"4l) apply to the full range of the variables, i.e. they 
represent the kernels completely. The function V qq t (t, t', r) was already derived in refs. |2D|, [13], 
the complete treatment of both singlet cases are given in |J . For the calculation one may fix the 
parameter r, apply the scaling properties and take into account the necessary general structure 
in the t, t'-plane derived in Note, that our kernels are represented in a very compact form. 
Equivalent results obtained by other authors [9-11] consist at least of two separately calculated 
expressions. For r = 1, our expressions correspond to £ 
our parameter r equals to — £/2. Note also refs. [12 



1 in 



In the notation of ref. [11 



6 Special Cases 

The evolution kernels given above cover limiting cases which were studied before. These are 
characterized by special kinematic conditions for the matrix elements. This concerns the case of 
forward scattering (p 2 | — ► {pi\ = {p\ and the transition from the vacuum state (0| to a hadron 
state {p\ wave functions. 

6.1 The Brodsky— Lepage Limit 

For r = ±1 the equations (|38|-[54]) correspond to the limit (pz\ —>■ (p\, (pi\ — > (0|, which is known 



as the Brodsky-Lepage |L4[ and Efremov-Radyushkin |L5[ case. This limit p\ — ► may be 
performed formally, i.e. irrespective of any other quantum numbers, leading to correct results, 
cf. [p0 |. The corresponding evolution equations are given by eq. (0) using as variables x and y, 



given in eq. (|3"T|). As an example we consider the simplest case in the range x, y < 1, x ^ y 

V«(x,y) = c F {Q(y -x)[---+ — Ur] + Q(x - y )[\^- -^—+ r^]}- (55) 

y V (y-x) 1-y 1-y (x-y) 

6.2 The Near Forward Representation 



This representation has been introduced in [|Tl|] and also used in ref. 0. Essentially it consists of 
the part for t > r and t' > r of the general kernel. It contains the forward scattering evolution 
kernels case as limiting cases. Especially, for a correct application of the evolution equation to 
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the near forward matrix elements this representation is needed (for t > r ) but additionally 
another representation for t < t has to be taken into account. In an example we show that this 
representation follows from our general kernel, i.e. that the general structure of V qq t covers also 
the case t > T,t' > r. For this range sign(y) = — sign(y) and 0(1 — |) = Q(y — x) holds. Using 
these changes and dropping, for simplicity, the +-prescriptions we obtain 

V»(z,y) = C F e(y-s){-[l + — ]-^— + — ]} 

y y — x 1 — y x — y 

= C F e(y-x) — [l + ^}, (56) 
y-x yy 



where again x and y are taken from eq. (|37j) 



6.3 The Altarelli— Parisi Limit 

The Altarelli-Parisi kernels can be directly determined from the anomalous dimensions and 
suitably defined quark and gluon distribution functions 

P(z,fi) = — / d(2px^)(p\O q \p)(^, f i) — — , (57) 
Ztx J-oo 2ipx 

zf G (z,fi) = —J^d{2p^){p\0 G \p){K.^)-^-^. (58) 

The respective polarized parton densities are obtained by replacing f q ' G (z) by Af q ' G (z) and O q ' G 
by Of G . In ref. this procedure has been described and leads to the known results, 0. On 
the other hand we should remark, that our general expressions for the one-variable evolution 
equations und the corresponding kernels contain the Altarelli-Parisi case in the limit r — ► 0. The 
determination of this limit is nontrivial because cancellations of divergent terms occur. As an 
example we formulate the basic steps for the kernel V^f. 

(2x-l)V G J(x,y) = CV{e(-)©(l--)sign(y)(l--) 

y y y 

-©(=)©(! -I)sign(y)(l- (59) 

y y y t 

We re-express the variables x, y again by the t, t'-variables. Finally the limit r — >• is performed 
yielding 

limK^(M') = 1,C F { 1 + {1 ~ Z)2 \ z = * (60) 

This is the corresponding Altarelli-Parisi splitting function. Another possibility consists in de- 
riving first the near-forward representation and performing then the limit to the forward case, 
cf. H. 



7 Solutions 

Here we look for solutions of the two variable equations in the singlet case. We extend the 
method of ref. [17| to this case. The idea follows the original solution of the Brodsky-Lepage 
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equation given by Efremov and Radyushkin [fL5 |. The method consists in forming moments of 
the partition functions such that the resulting evolution equations for these moments can be 
solved. We start form eq. (p7]) and look for equations for the moments of the partition functions 

(0), 



Gl(z. 



dz+z+ G (^_|_, z_ 



(61) 



One obtains from eq. fl27j) 



d r+i 



dfi 2 J-i 



dz^.z^. G (^+ ? z_ 



2tt 2 7-i 



+1 



- 1 dz\ 



1 z+ Vz+ 



K lJ \ Z-- ^yz'-,^y \ G ij (z' + , z'_) (62) 

\ Za_ Za_ I 



la s (/i 2 ) f +1 dz + ( z + 



n—hi 



2 2tt 7-1 Izi 



+1 



dz' 



+1 



yielding 



with 



1 V 4 4/ 



2 2tt 



2\ ,+1 



-1 



(64) 



T^(z.,z'_ 



-1 V^l 



n—hi 



K^{z— Z -±z>J-±) 

Z_L Z, 



(65) 



In this way a first diagonalization is obtained. The kernels r^(z_, 2/.) for the unpolarized case 
read 



c i? {e(z-,')[(i-^r(i+ 2 
+e(/-z)[(i^r(- + 



z — z' 
2 



I + 2' n (2;' — z) 



)] + &(*-*')}, 



-2N f T R {e(z - m\^rX- + 2,1(1 " z ' ] - {2zz ' - 



[n + l)n(n — 1) 



+e(z'-z)[( 



l + z n 2 + 2n(z' - z) - (2zz' - 1) 



]}• 



(66) 



l + z ,J (n+l)n(n-l) 

-C F {Q{z-z')[{\^ 

1 — z' n 

+e(z'-z)[(^r^+S(z-z')}, 

C A{ e(z - , 0[( l^ r( 3(l + + g^3(l + z Q(l -z) _ 2 + 



(67) 



(6? 



1 - z" 



n — 1 



n + 1 



AM- U ,1 + 1 W 3(1 Z ^)(1+^) , 6*^3(1 -^(l + z) 2 
WlZ jK l + ^ j 1 n-1 + n ~ 



n [z — z 
2 



n + 1 



n (2;' — z)_ 



]} 
(69) 
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Similarly one obtains for the polarized kernels 



C F {Q(z-z')[( 



1 - z 
1 - z' 



- + 7- 

*n (2 — 2 



AT G q{ 



z, z 



AT° G (z,z' 



+Q ( Z ' _ z )[(i±£)»(I + 2 

V nK l + z fJ K n (z'-z) + 
-2N f T R {e(z - z')[^rl 

+C F {Q(z - z f )[(^r-\ 
1 — z' n 

c A {e(z - z ')[{l^n- + 



)] + 38(z-z')}, 



1 - z" 



1 + z 4 2 

e(2'-2)[(-^4r(-+— — 

v m l + z 1 ' K n (z' - z 



a ( : — z') + 

)}} + f3 Q S(z-z'). 



(70) 



(71) 



(72) 



(73) 



As a next step we have to find a solution of eq. (|B5|). One finds that the following symmetry- 
relations hold for kernels T GG (z, z') and AT GG (z, z') : 



(i- z > 2 rri?(z,z') 

'l-z >2 ) n ATtiHz,z') 



r^( z ',z)(i- z *r, 

AT^(z',z)(l-z 2 r 



(74) 
(75) 



This is expected because of similar problems studied previously |17|, [T5J. The final diagonalization 
of eq. fl64|) , referring to the partition functions eq. (|l0j) for which the quarkonic and gluonic 
operators are dealt with equally, can thus be performed using Gegenbauer polynomials. 



- E rK»(i - z 2 rcz +1/2 (z)cz +1/2 (z')N n>m , 

The coefficents N n ^ m are the normalization factors of the Gegenbauer polynomials 



(76) 
(77) 



+i 



dz(l - z 2 rC? +1/2 (z)Cr ir \z) = S lk N n , k = S, 



n+1/2. 



Ik 



)2n+l r 2 



T z n + 



1\ (k + n + 1/2) • k\ 



i -i 



2 J 2irT(2n + l + k) 



in 



r^ J m and Ar^ m are the respective expansion coefficients which can be easily calculated. As 
well-known [^TJ the various possible representations of the solution of evolution equations using 
expansions in orthogonal polynomials are badly convergent p2|, E3J requesting high-precision 



representations |23| which are difficult to handle in practical applications. 

Another method of solution consists in a forming also Mellin moments in the variable Z-, 
f+L dz-ztG l n (z-) = G l nk . Since the kernels T)?(z-, z'_) and AY${z-, z'_) obey the representation 



r n i(z_,z'_)= I dw 2 w n 2 - h3 K i ' j ^--w 2 z'^w 2 ). 



the evolution equations may be written as (see also 

a s (/i 2 ^ 



2 d r 



2 nk 



2tt 



1 n,kl nl ■ 



(79) 



(80) 



3,1 
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with the transformed kernels given by 

= U dw wT l - hj w\- l k^{w u (81) 

These evolution equations are not diagonal with respect to the indices (k, I). The explicit ex- 
pressions for T l ^ kl and are given in ref. [7,5b]. For fixed n they form triangular matrices. 
The eigenvalues are the diagonal elements k = I 

Kb* = \ J Dw w n 2 +k - hj K*' j (wi, w 2 ) = ^ k _ hj . (82) 

The coefficients r i^ +k _ hj are the anomalous dimensions of the forward case with a shifted Mellin 
index. 

8 Conclusions 

The evolution kernels for the twist 2 light-ray operators both for the case of unpolarized and 
polarized deep inelastic non-forward scattering were derived for the flavor non-singlet and singlet 
cases. In general the partition functions depend on two distribution variables. One may study 
as well specialized evolution equations in one distribution variable implying external constraints, 
covering the case of evolution equations for non-forward parton densities. In this way, among 
various others, also the well-known evolution equations as the Brodsky-Lepage or Altarelli- 
Parisi equations can be obtained. The solution of the two-variable evolution equations in the 
non-singlet and singlet cases can be either performed applying a single Mellin transform and 
using Gegenbauer polynomials or by a two-fold Mellin transform. 
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